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ABSTRACT 

GRB 080607 has provided the first strong observational signatures of molecular absorption bands to- 
ward any galaxy hosting a gamma-ray burst. Despite the identification of dozens of features as belong- 
ing to various atomic and molecular (H2 and CO) carriers, many more absorption features remained 
unidentified. Here we report on a search among these features for absorption from vibrationally- 
excited H2, a species that was predicted to be produced by the UV flash of a GRB impinging on 
a molecular cloud. Following a detailed comparison between our spectroscopy and static, as well as 
dynamic, models of H2* absorption, we conclude that a column density of cm~^ of H2* 

was produced along the line of sight toward GRB 080607. Depending on the assumed amount of 
dust extinction between the molecular cloud and the GRB, the model distance between the two is 
found to be in the range 230-940 pc. Such a range is consistent with a conservative lower limit of 
100 pc estimated from the presence of Mg I in the same data. These distances show that substantial 
molecular material is found within hundreds of pc from GRB 080607, part of the distribution of clouds 
within the GRB host galaxy. 

Subject headings: galaxies: ISM — gamma rays: bursts — ISM: clouds — ISM: molecules — molecular 
processes 



1. INTRODUCTION 

Long-duration GRBs are currently understood to be 
the electromagnetic manifestation of highly beamed en- 
ergy from the sites of the demise of young massive 
stars (see, e.g.. IWooslev fc BloomI [20061 ) . This implies 
that GRB progenitors must have formed and lived near 
molecular clouds in their respective host galaxies. Di- 
rect evidence for the molecular nature of their birth- 
place has been lacking, presumably owing to complete 
destruction of molecules within ~100 pc o f the GRB 



dTumUns on et all 120071: IWhalen et all l2008f l . Recently, 
IProchaska et al.l ( 20091 hereafter P09) presented the 



first unambiguous spectroscopic detection of absorption 
bands from H2 and CO molecules along a translucent 
GRB sight line (GRB 080607, z = 3.0363). Their ini- 
tial analysis provided the first discovery of molecular gas 
with properties very similar to those of Galactic molecu- 
lar clouds, but located in the star-forming ISM of a GRB 
host galaxy. 

The detection by P09 of neutral species (e.g., C I, Mg I) 
toward GRB 080607 indicates that the bulk of the atomic 
gas was >100 pc from the GRB, presumably beyond 
the photoionization sphere formed by t he progenitor 
( Whalen et al.l [20081 ). However, whereas IWhalen et al.l 
()2008f l introduced an additional galactic-wide FUV field 
as means to suppress the formation of unobserved H2, 
the discovery of more than 10^^ cm~^ of H2 toward 
GRB 080607, well shielded by dust extinction (rest-frame 
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Av = 3.2 mag), raises the possibility that the molecular 
material could have been in close proximity to the GRB. 
The detection by P09 of absorption from rotationally- 
excited CO, and the determination that atomic and 
molecular lines are kinematically distinguished by 30 ± 

1 5 km s~\ spec ifically allow for this scenari o. 

iDraind (|2000D . followed bv lDraine fc Had (|2002l . here- 
after DH02), predicted that observable quantities of 
vibrationally-excited H2 (from v" > 0, hereafter H2*) 
can be produced by UV photons in molecular gas within 
100 pc from a GRB source. With an assumed model 
peak luminosity of Lq = 1/40 of that observed for GRB 
990123, DH02 produced a column density [N) of wlO^^ 



of H2* 



in a gas with an initial N — 10 cm 



of cold H2 at close proximity (<1 pc) to the GRB. P09 
reported cold H2 with N « lO^^-^ cm'^ toward GRB 
080607, suggesting the possibility that H2* could also 
have been present. 

Here we report on our successful search for absorption 
from H2* toward GRB 080607. Previously, one Galactic 
sight line was sho wn to have ^500 a bsorption lines asso- 
ciated with H2* (jMever et al.ll2001l . toward HD 37903), 
while the sample of low-redshift galaxies, meanwhile, 
that are found to emit ro- vibrational H2* lines near 2.1 
/im has been inc reasing steadily , albeit m ostly attributed 
to shocked gas ([Th ompson et al. 1978; J affe fc Bremeil 
119971 : iDonahue etal.. 2000; Dale ct al. 200"9[). Our sight 
hue toward GRB 080607 at the redshift of z = 3.0363 is 
thus the most distant detection of H2*, as well as the first 
extra-galactic detection of H2* via absorption: this was 
made possible thanks to an extremely bright GRB, an 
extremely capable optical system (the LRIS spectrome- 
ter on the Keck I telescope), and a rapid response to the 
GRB alert from Swift. We also show, by a comparison 
with theoretical predictions, that this detection serves to 
confirm the general picture of the interaction between a 
GRB and its immediate galactic environment. 
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Fig. 1. — Zoomed comparison between the Keck/LRIS R1200 data (black histogram) and the H2* dynamic model (smooth red line). 
Brown masking designates absorption from other identified species (e.g., Si II at 6163 and 6189 A), where the detection of H2* would be 
compromised by blending. The agreement is very good between stronger H2* absorption features and equally strong, hitherto unidentified 
"lines" in the spectrum. The noise level is denoted by the green line at the bottom. 



2. OBSERVATIONS AND MODELING 

The acquisition of the spectroscopic data analyzed here 
was described by P09. Here we recap some essential in- 
formation relating to the sequence of events in the rest 
frame of the GRB, owing to the predicted temporal evo- 
lution of the H2* abundance (DH02). 

GRB 080607 trigg ered the Swift detect or at 06:07:27 
UT on 2008 June 7 (|Mangano et al.ll2008[ ). here defined 
as i = s. A first series of three Keck/LRIS exposures 
was taken by us (DAP and JSB) with the B600 and 
R400 gratings, simultaneously employing the blue and 
red cameras. These exposures were 75, 150, and 300 co- 
moving seconds long, centered on rest-frame times t = 
340, 490, and 730 s. A second series of two exposures 
was taken with the B600 and R1200 gratings, both 370 
s long, and centered on the post-trigger times t = 1230 
and 1640 s in the rest frame of the GRB. In terms of 
spectral resolution, the sub-arcsecond seeing resulted in 
A/AA = i? - 4000, 2000, and 1000, for gratings R1200, 
B600, and R400, respectively. 

Our initial search for possible H2* features involved 
spectral syntheses of relevant transitions with the Y.S. 
code Ismod.f, and comparisons with the highest-i? data 



from the R1200 grating. Each transition was modeled 
with a single Voigt profile, employing as fixed parameters 
the resolution R, transition rest wavelength (Arcst) and 
oscillator strength (/-value), the Doppler parameter b, 
and a radial velocity of +30 km s"""^ relative to the atomic 
gas (P09). Both the Arcst and /-values for B — X and C — 
X transitions of H2* were downloaded from the MOLAT 
website® (E. Roueff, personal communication). Table 1 
lists our H2* identifications for 25 absorption features 
detected at the >4-tT level in the R1200 spectrum. 

Once we were convinced of the presence of (static) H2* 
absorption in the data toward GRB 080607, we shifted 
to full dynamic modeling, with photoexcitation, pho- 
todissociation, and photoionization of the gas treated 
following DH02. The UV and X-ray emission from GRB 
080607 was assumed to be described by 



13.6 eV 



1+/3 



4{t/to 



(1) 



[1 + (t/to 

where Lq = 5 x 10^*^ erg s~^, /? = —0.5, and tQ = 10 s. 
With Ay = 3.2 mag (P09), this reproduces the observed 
K-band fiux density at iobs = 300 s. The value of Lq is 

^ http://molat.obspm.fr/ 
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TABLE 1 

Features Identified toward GRB 080607 





Wx'' ID in P09'' 


-^rcst 


ID in DH02'^ New ID'^ 


Nonni 


(A) 


(A) 


(A) 


(f"; J") 





5762.6 
5783.5 
5789.7 
5797.4 
5988.6 
6002.3 
6012.5 
6115.3 
6124.8 
6130.3 
6135.6 
6156.6 
6179.2 
6205.8 
6213.8 
6263.7 
6269.1 
6318.8 
6414.3 
6432.1 
6477.4 
6489.1 
6525.8 
6604.4 
6614.6 



0.8 
1.1 
0.6 
1.4 
1.0 
2.6 
1.9 
1.1 
1.5 
1.3 
0.9 
0.9 
0.6 
0.6 
0.6 
1.4 
0.6 
0.7 
0.9 
0.7 
1.6 
1.5 
0.6 
0.7 
1.6 



UID 
UID 
UID 



w/ S 



II 



w/ CIV 

UID 

UID 
w/ Si I 

UID 

UID 
w/ Fe II 

w/ Fe II 
w/ Fe II 



1427.1 
1432.3 
1433.8 
1435.7 
1483.1 
1486.5 
1489.0 
1514.5 
1516.8 
1518.2 
1519.5 
1524.7 
1530.3 
1536.9 
1538.9 
1551.2 
1552.6 
1564.9 
1588.6 
1593.0 
1604.2 
1607.1 
1616.2 
1635.7 
1638.2 



Pl+Rl 
RO+Rl 



Rl+Pl 
PI 



RO+Rl 
P1+R3 
R1+R2 



(6; 0,1) 
(6; 2) 
(6; 4) 
(8; 0) 
(10; 2) 
(7; 0,1) 
(7; 1,2) 
(8; 0,1) 
(7; 2) 
(11; 0,1) 
(8; 2) 
(7; 3) 
(10; 1) 
(10; 3) 
(9; 0,1) 
(12; 1,2) 
(13; 1) 
(13; 1) 
(10; 0,1) 
(10; 2) 
(9; 1) 
(9; 2) 
(12; 5) 

(9; 1) 
(11; 2) 



Detections with Wx > 0.6 A (>4(t).'^ UID: Listed in P09 
without ID. Otherwise, was included within a stronger atomic 
feature.*^ Strongest transitions within 1.0 A of Arost for R = 
350. Strongest contributing ro-vibrational levels for R = 4000. 

200 times higher than that adopted in DH02. Indeed, 
GRB 080607 ranks among the most p owerful to date, a n 
equal to GRB 990123 in 7 luminosity (|Golenetskiill2008D . 

Total line of sight abundances for H and H2 were taken 
from P09, with the H nucleon density (nn) taken to 
be 10^ cm~^. The dust/gas ratio in the pre-GRB gas 
was assumed to be ~10% of the local Milky Way value 
(P09). For photoexcitation of the H2 and photoioniza- 
tion of H and H2, we were primarily interested in the 
dust extinction cross section per H nucleon at A < 1110 
A, which we take to be cr = 2 x lO"'^^ cm~^, so that 
A ' = 5.33 X 10^2 X 2 X IQ-^a /i 086 = 9.8 mag, or 

A '/Av^S.l. 

lOOOA' 

Finally, the evolution of the irradiated dust was fol- 
lowed with the same assumptions as in DH02. Initial 
grain size was a = 0.04 /im, which could decrease by 
thermal sublimation if the grains became hot enough. 
For calculating UV extinction we used Qext = 2, where 
Qoxt is the ratio of the extinction cross section to the 
geometric cross section, 7ra^. 

3. RESULTS AND DISCUSSION 

3.1. H2 * Column Density 

Fig. [H shows a small portion of the R1200 data, super- 
posed by 336 modeled transitions of H2*. This region 
includes five absorption "lines" that were listed, but not 
identified, in P09, four of them as a group between Aobs 
= 6115-6130 A. However, each H2* "hue" is a blended 
feature that includes several strong transitions. It is ev- 
ident that the four observed features near 6120 A are 
too strong to be attributed to noise: their equivalent 
width {Wx) values listed in P09 range over 1.1-2.3 A, or 




6600 6650 
Wavelength (Ang) 



Fig. 2.— Same as Fig. [T] but showing 1000 A of the R1200 
data and the model of H2*. Dozens of narrow H2* and observed 
features appear to agree with each other. The global pattern of 
variable v' — v" band strength is also matched. 

7-15 times the l-cr detection limit. Since these "lines" 
cannot be identified with any other plausible atomic or 
molecular species, at the redshift of either GRB 080607 
or the two intervening absorption systems, we are con- 
fident that H2* provides the only viable model for their 
presence in our data. For the best match with the data, 
the static modeling with Ismod.f returned log iV(H2*) 
= 17.7 ± 0.3, using b — 2 km s^^, a value taken from 
the CO analysis of P09. Our full dynamic modeling with 
the DH02 code, computing H2 excitation for 10'^ s in the 
GRB time frame and using 6 = 3 km s~^, returned log 
iV(H2*) = 17.5 ± 0.2, in good agreement with the static 
model. This served to confirm that the first case of H2* 
excitation by GRB photon pumping has been identified. 

The entire range of the R1200 spectrum (Arcst = 1400- 
1720 A) covers the red portion of H2* absorption, and 
includes 3372 transitions with J" = 0-29, v" — 1-14, and 
v' = 0-16. Figure [2] presents a global view of the H2* 
model spectrum superposed over 1000 A of the R1200 
data. The overall agreement between data and model is 
acceptable, given that the S/N « 10 and the presence 
of numerous other absorption features from, e.g., host- 
galaxy CO bandhcads and Fe II from two foreground 
absorbers (P09). The few instances where the modeled 
absorption is significantly deeper than the data involve 
levels with v" > 9, presumably related to our incom- 
plete description (underestimation) of UV photoioniza- 
tion rates out of high-w" levels. 

We identified 25 narrow H2* features that are the least 
blended with other species and list these in Table [T] Pre- 
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Fig. 3. — Same as Fig. (2] but showing the lower-resolution R400 
data and H2* model. Absorption from wide features made up of 
numerous H2* transitions is determining the shape of the contin- 
uum of the R400 spectrum. Some absorption features previously 
assigned to atomic carriers are clearly significantly blended with 
contribution from H2*. 

viously, P09 estimateci atomic abundances in the weak 
limit from measure(i W\ of absorption features. Since 
some transitions listed in P09 are seen here (Fig. [2|) to 
be blended with some of the stronger H2* features, cer- 
tain Wx values should be revised downward. For exam- 
ple, the H2* feature at 6466 A falls into a blend with the 
Ge II 1602 line, contributing 2.0 A to the total Wx of 5.3 
A and reducing its weak-limit abundance (estimated by 
P09 to be super-solar). However, only a future detailed 
study of sight-line saturation will be able to provide ro- 
bust abundance determinations based on H2*-corrected 
data. 

The presence of H2* in our R1200 spectrum is con- 
firmed by the lower-S/N exposures from the B600 grat- 
ing {R — 2000). In all, the blue and red gratings provide 
~500 A of rest-frame coverage of numerous absorption 
features belonging to (or blended with) H2 * toward GRB 
080607. 

3.2. H2 * Variability vs. Non- Variability 

The theoretical work of DH02 followed the temporal 
behavior of H2* abundance in response to the initial UV 
flash from a GRB afterglow. Upon photoexcitation into 
the electronic B and C states, only ~15% of H2 is pho- 
todissociated, with the majority, therefore, decaying im- 
mediately into w > levels of the X ground state. A 
very rapid rise in A^(H2*) occurs over the time scale of 
the burst («10 s) as the UV photons get absorb ed by 
the cold H2 gas. The 10^ s lifetime of H2* (|Drainell200a) 



means that it remains around for almost two months in 
the observer frame. 

As described in §2, our two groups of red exposures 
were centered on i ~ 500 and 1400 s in the GRB rest 
frame. As shown in Fig. [3l a degraded R — 1000 ver- 
sion of the H2* model that matches the R1200 spectrum 
shows remarkable agreement with the upper envelope of 
the R400 spectrum. Thus there is no indication of varia- 
tion in A^(H2*) over the At ~ 900 s that elapsed between 
the clusters of R400 and R1200 exposures. This obser- 
vational result confirms the DH02 theoretical predictions 
that iV(H2*) varies over time scales that are either signif- 
icantly shorter or significantly longer than those covered 
by our observations. 

3.3. The Distance of H2 * from the GRB 

The value of A^(H2*) detected here is only ~3% of that 
produced by the models in DH02. According to DH02, 
the amount of H2* produced depends on the fluence, F, 
of the source, and hence on the ratio Lo/D^, where D is 
the distance from the GRB to the center of the initial H2 
integration zone. Owing to the much higher value of Lq 
used in our modeling of GRB 080607 (§2), the production 
of an appreciably sma//er amount of H2* relative to DH02 
would require a very large increase in D. 

We employ the detailed modeling of DH02 in order to 
explore the dependence of A^(H2*) on D. The effect of 
the GRB on the H2 will depend on what fraction of the 
H I is located between the GRB and the H2 , because the 
dust in the H I ^as can help shield the H2 from photons 
with A < 1110 A, and the H I itself will help protect the 
H2 from photoionization. 

We consider two cases. Model A has iV(H°) = 2 x 10^^ 
cm^^ located between the GRB and the H2, with the 
remainder of the H I located beyond the H2. Model B 
has N(H°) = 5 X 10^2 cm'^ (i.e., all of the H I from 
P09) located between the GRB and the H2. For a given 
large value of D, model B will have less UV-pumped H2 
than model A, because of the additional dust shielding 
in model B. However, for smaller values of D, reduced 
shielding by the destruction of dust and H I in model 
B may result in complete destruction of the H2 by pho- 
todissociation or photoionization. 

The upper panel of Fig. |4] shows iV(H2*) versus D for 
the two models. The observed value of w 10^^'^ cm^^ 
(§3.1) is generated for D = 2.9 x lO^^ or 7.3 x lO^" cm 
for Model A or B, respectively. Thus we conclude that 
the H2 in GRB 080607 is located between 230-940 pc 
from GRB 080607, a value much larger than that used 
in DH02 {D < 1 pc). 

The lower panel of Fig. [4] shows the fraction of the ini- 
tial dust and H2, that survives after the GRB flash. For 
more than 50% of the initial grain mass to be sublimed, 
the H2 gas would have to be within ~190 pc, slightly 
closer than the location of the H2* in model B, and sub- 
stantially closer than in model A. We conclude that for 
these two models, the bulk of the dust survives the GRB 
flash, consistent with Ay « 3.2 mag observed during the 
afterglow phase. Substantial destruction of the H2 by 
a combination of photodissociation and photoionization 
occurs only for D < 2 x 10^° cm = 65 pc. Given the ob- 
served A^(H2*)/A^(H2) « lO^^ "^, we conclude that there 
has been minimal destruction of H2 in the cloud where 
the H2* resides. 
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Fig. 4. — Column density of vibrationally-excited H2 as a func- 
tion of distance from GRB 080607, for two models differing in the 
location of the H I relative to the H2 (see text). 

More complex models are of course possible: if more 
than one H2 cloud is present, there could have near- 
complete destruction of the H2 in the inner cloud, leaving 
only a small amount of H2 with high levels of vibrational 
excitation, with the bulk of the H2 in a more distant 
cloud with little vibrational excitation. 

3.4. What about CO*? 

The presence of H2* prompted us to search for 
vibrationally-excited CO (CO*) in our data. P09 re- 
ported the detection of log N = 16.5 ± 0.3 cm~^ of cold 
(f " = 0) CO, all of which was rotationally-excited (here- 
after r-CO*) at Tex 300 K up to J" = 25. If the vi- 
brational excitation of CO owing to UV pumping scales 
with that of H2, then we should expect log A^(CO*) — log 
A^(CO) - 3.7 « 12.8 cm-2. CO* transitions for v" > 
bands were taken from Kurucz line lists^ and modeled 
with Ismod.f, but no match with the R1200 data was 
found. Based on the strongest v" = 1 band at 6400 A, 
a feature with A^(CO*) = 1 x 10^^ cm~^ would have a 
depth of 10%, comparable to the 1-ct noise level in the 
data. 

We note that the photophysical behavior of CO dif- 
fers from that of H2 in two important respects. First, 
owing to the existence of strongly predi ssociating Ryd- 
berg states in CO below Aicst = 1100 A (jLetzelter et al.l 

http://kurucz.harvard.edu/LINELISTS/LINESMOL/ 



[19871 Ivan Dishoek fc Blackl [Tosi: iSheffer et all [200l . 
CO*/CO should be much lower than H2*/H2. Second, 
CO* spontaneously dec ays via dipole tran sitions with a 
hfetime of ca. lO'^ s (lO kada et ani2002[ ). roughly 10^ 
times faster than the quadrupole transitions of H2*. Any 
small amount of surviving CO* would be rapidly con- 
verted into r-CO*. However, while CO photophysics may 
explain the absence of CO* toward GRB 080607, the 
challenge of accounting for the conversion of the bulk of 
cold CO into r-CO* remains. 

Both dust extinction and cold H2 provide very effective 
UV shielding of CO below Arest = HOO A, allowing it to 
survive in deeper layers of a cloud. The non-dissociating 
A~ X bands (with v' < 8) are significantly less shielded: 
the data toward GRB 080706 (Fig. 1 of P09) show that 
the UV flux level at Arost > 1322 A is ~5 times higher 
than at Arcst < 1150 A. UV pumping into the CO A state 
could result in its complete conversion into CO* and then 
into r-CO*. However, assuming that initially the CO 
is cold (Tox 11 K, J" < 5), the production of r-CO* 
that happens to maintain a quasi-thermal distribution at 
^300 K up to J" — 25 must be considered a coincidence 
of indirect excitation. 

Radio observations show that collisionally- 
excited r-CO* can be fo und in Galactic PDRs 
(iTielens fc HoUenbachl with T^^ ^ 100-1000 

K (Ivan Dishoek fc B lack T98| iDraine fc Bertoldil [l996l : 
iHollenbach fc Tielen£i,1997j) . This range agrees with Tex 
~ 300 K for r-CO* toward GRB 080607, subject to the 
caveat that appreciably higher values of A^(CO) and Ay 
are found in Galactic PDRs. We surmise that the r-CO* 
toward GRB 080607 could be a pre-burst observational 
signature of the PDR produced by the progenitor of 
the G RB, as it was formin g a giant H II region around 
itself (|Whalen et all l2008f) . A detailed calculation of 
this process is clearly warranted. 

4. CONCLUDING REMARKS 

In addition to providing the first positive detection 
of H2 bands and the first observation of CO bands in 
a GRB- host galaxy (P09), we have shown that GRB 
080607 also provides the first evidence for H2* in a GRB 
host galaxy and marks the highest-redshiftcd H2* de- 
tected to date. 

This discovery of UV-pumpcd H2* toward GRB 080607 
serves to confirm the predictions for its produc tion un- 
der precisely such circumstances ([Draind I2OOOI DH02). 
Our initial static modeling with Ismod.f indicated log 
iV(H2*) < 18.0 toward GRB 080607, significantly lower 
than the original predictions of DH02. We then showed 
that the dynamic models of DH02 can successfully repro- 
duce observed H2* absorption once modifications involv- 
ing individual GRB luminosities and adjustable distance 
to the molecular cloud are incorporated. Thus the DH02 
photoexcitation code shows that the bulk of the molec- 
ular cloud harboring H2* is located at a model distance 
of 230-940 pc along the line of sight from the UV after- 
glow. This is much farther than the original arrangement 
in DH02, but still in the local galactic neighb orhood of 
GRB 080607. In the scheme of IWhalen etall (|2008f ). a 
model distance of >230 pc for the molecular cloud means 
that it is located outside the «100 pc radius of H II re- 
gion carved by the progenitor of GRB 080607. 

One interesting result of the increased distance from 
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the GRB to the gas is the inabihty of the radiation beam 
to destroy dust embedded in the cloud. Whereas in DH02 
dust was easily destroyed at close quarters to the GRB, 
following its heating to 2000-3000 K, this process can- 
not operate at larger distances, leading to survival of the 
dust in the cloud. Such dust remains on the line of sight 
between the GRB and the H2 destruction front, dimin- 
ishing the effects of such destruction. 
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